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Abstract 

Using a recently developed bootstrapping method, we compute the phase dia¬ 
gram of scalar field theory on the fuzzy disc with quartic even potential. We 
find three distinct phases with second and third order phase transitions between 
them. In particular, we find that the second order phase transition happens 
approximately at a fixed ratio of the two coupling constants defining the poten¬ 
tial. We compute this ratio analytically in the limit of large coupling constants. 
Our results qualitatively agree with previously obtained numerical results. 


1. Introduction 


By a fuzzy space, one usually means a geometric quantization of a compact Kahler man¬ 
ifold. The compactness of the space implies that the arising Hilbert space and therefore 
also the algebra of observables given by the endomorphisms of this Hilbert space are fi¬ 
nite dimensional. Correspondingly, there is a minimal resolution with which the space is 
perceived, which renders it “fuzzy.” For this reason, fuzzy spaces are good candidates for 
regularizing quantum field theories [1], because the path integral over all observables turns 
into a finite number of ordinary integrals. To compare to what extent quantum field theo¬ 
ries on fuzzy spaces provide an approximation to the corresponding quantum field theories, 
it is particularly useful to study the phase diagrams of the theories in the thermodynamic 
limit. 

To this end, we need to evaluate the free energy of the fuzzy quantum field theory. As 
usual in geometric quantization, real functions on compact Kahler manifolds are mapped 
to hermitian operators, representing the quantization on the corresponding fuzzy mani¬ 
fold. Therefore, scalar field theories on fuzzy spaces are simply hermitian matrix models. 
Contrary to the hermitian matrix models most common in the literature, however, these 
matrix models come with a kinetic term containing fixed external matrices. This kinetic 
term presents an obstacle for applying the usual techniques for solving matrix models di¬ 
rectly. In particular, the action is no longer invariant under similarity transformations and 
diagonalizing the matrix is no longer straightforward. This problem can be overcome by 
rewriting the kinetic term as a multitrace expression, which can then be solved at least in 
the limit of large matrix size, e.g. by the saddle point approximation. 

This approach was developed in [2] and used there, in [3] and in [4] to compute the 
phase diagram of scalar quantum held theory on fuzzy complex projective spaces as well 
as on a three-dimensional space consisting of the Cartesian product of R and the fuzzy 
sphere. Further applications of this technique are found in [5], see also [6, 7, 8]. The 
rewriting of the kinetic term was done by applying techniques from group theory, making 
the calculations somewhat cumbersome. An alternative bootstrapping method for turning 
the kinetic term into multitraces was then found in [9]. Here, enough conditions on the 
coefficients in the multitrace expressions are derived to fix them uniquely. 

The purpose of this letter is to use the bootstrapping approach to compute the phase 
diagram of scalar field theory on the fuzzy disc [10] and to compare the result to the 
numerical findings of [11]. The fuzzy disc is particularly appealing as the kinetic term is 
somewhat simpler than in the case of the fuzzy sphere. One may therefore hope that on 
the fuzzy disc, quantum scalar field theory is better behaved than on the fuzzy sphere or 
even that the resulting hermitian matrix model is fully solvable. 

This letter is structured as follows. In section 2, we briefly review the construction of 
the fuzzy disc and the definition of quantum scalar field theory on it, fixing our conventions. 
Section 3 deals with rewriting the resulting action as a multitrace expression and taking the 
limit of large matrix size. We compute the phase diagram of the model using a saddle point 
approximation in section 4, where we also compare our results to the numerical literature. 
We conclude in section 5. 
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2. The model 


We start our discussion with a concise review of scalar field theory on the fuzzy disc. 


2.1. Fuzzy disc 

The fuzzy disc [10] provides a quantization of the algebra of functions on the unit disc in 
C = R 2 . It is obtained by truncating and rescaling the matrix algebra of the Moyal plane. 
Moreover, one can obtain the fuzzy sphere by gluing together two fuzzy discs. Below, we 
briefly recall its construction, following roughly [10]. 

Recall that the Moyal plane Rg is the geometrical quantization of the complex plane 
with respect to the canonical symplectic structure, see e.g. [12] for details. The result 
of this procedure is an infinite dimensional Hilbert space H which agrees with the usual 
Hilbert space of the harmonic oscillator up to a normalization. That is, we have a vacuum 
state |0) together with annihilation and creation operators satisfying 

a|0) = 0 , [a, at] = 9 , 9 € R >0 . (2.1) 

We denote the eigenstates of the number operator N = ad) with eigenvalue n6 by |n): 

!V|n) = n0\n) , n € N . (2-2) 


We also introduce a corresponding basis for P*: (m\, m € N, normalized such that (m\n) = 
5 mn . The endomorphisms on H are spanned by linear combinations of the operators |m)(n|. 
The Berezin symbol map a arising in the quantization assigns to each element of End('R) 
a function on <D = R 2 , providing a dequantization map. We have 


OO 

f = fmn\m)(n\ , 

m,n =0 


v(f)(z,z) 



°° ~m „n 

V f __ 

^ mn V m\n\0 m+n 


(2.3) 


where z € (D. 1 In particular, as the usual quantization axioms demand, ct( 1) = 1. Note 
that real functions correspond to the hermitian endomorphisms End h(H)- 
Let us now introduce the following projector on a sub-Hilbert space H 0 . 


N—l 

P N ■= J2 \n)(n\ , U 0 := PnH and End(^ 0 ) := P N End(iR) P N . (2.4) 

n =0 

The Berezin symbol of this projector reads as 


a(P N )(z,z ) 


N -1 


£ 

n =0 


- e~~ 

n\9 n 


r(l + N/4) 

T(l + N) ’ 


r = \fzz . 


Here, T(n,x) is the incomplete gamma function. Note that 


Jim (a(P N )(z,z)) 

N^-oo 


1 for ^ < N , 
0 for ^ > N , 


(2.5) 


(2.6) 


i 


We follow physics conventions and write f(z, z) for a non-holomorphic function /. 
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which shows that the projector Pjv corresponds to a step function with support on the disc 
Dr C C of radius R = VN9. This justifies the identification of the algebra End#(P 0 ) 
with a quantization of the algebra of real functions on Dr, and we call the corresponding 
noncommutative space the fuzzy disc. We will always work with a disc of radius R = 1, 
fixing 9 = jj. 

2.2. Scalar field theory on the fuzzy disc 

To study scalar field theory on the fuzzy disc, we have to introduce a Laplace operator, 
i.e. an operator on End^("H 0 ), which approximates the usual Laplace operator A = 4^ J=. 
Recall that geometric quantization and the Berezin symbol map lead to the identification 

r\ r\ 

N[a, -] ~ gz and - N[af, ~ — (2.7) 

on the Moyal plane. On the fuzzy disc, we can combine these operators with the projectors 
Pjv to obtain a Laplace operator. There are two obvious candidates: 

A N f := -m 2 P N [a, [fit, f]\P N and A N f := -4 N 2 P N [a, P N [a f , f]P N ]P N (2.8) 

for / e End//(P 0 )- The first one was used in [10] and [11], but the second one has the 
advantage that AtvI = 0, as one might expect for the constant function a(l)(z,z) = 1. 
The latter expectation is somewhat debatable as constant functions on the disc are in fact 
step functions, and one might argue that the fuzzy boundary should lead to deviations from 
At\t 1 = 0. In the following, we will nevertheless work with the second Laplace operator 
for two reasons. First, this choice simplifies our computations dramatically and second, we 
will be mostly interested in the large N limit, in which both choices agree anyway. 

The second ingredient we need is the notion of an integral. Geometric quantization and 
normalization of the integrals lead to the following identification: 

J “ 2 “ = 7 ~§~ tr (/) = nd tr (/) = tr (/) > ( 2 - 9 ) 

where the last equality is again due to our choice R = y/N9 = 1. 

Introducing the shorthand a := PnclPn, we can now write down the action for scalar 
held theory on the fuzzy disc with quartic potential terms: 

S[« h] = ^ tr (—4 N 2 $[a, [a f , $]] + rl> 2 + g¥^j , (2.10) 

where r, g £ R such that rx 2 + gx 4 is bounded from below for all x € R. From now on, we 
will regard <3? € End h(Ro) as a hermitian N x IV-matrix 4> and drop the hat to simplify 
our notation. The partition function for the action (2.10) is then given by 

-2 := J dg D ($) e"^ , (2.11) 

where d/i/r( < l > ) denotes the usual Dyson measure on the space of hermitian matrices. 
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3. Computing the partition function 

As in the cases of the fuzzy sphere and fuzzy C P n in general [2, 3, 9], one cannot apply 
the usual techniques of hermitian matrix models to the partition function (2.11) in a 
straightforward manner. This is due to the fact that the kinetic term presents an obstacle 
to a simple diagonalization of <b. To circumvent this problem, we will rewrite the kinetic 
term as multitrace expressions. 


3.1. Multitrace action 


As usual in dealing with hermitian matrix models, we wish to diagonalize the hermitian 
matrix <J> as = flAfll, where A is a diagonal matrix containing the eigenvalues Ai,..., Xn 
of <J) and f l € U(./V). Under this change of variables, the Dyson measure d^£>(^ > ) on 
the space of hermitian N x A-matrices, which appeared in the partition function (2.11), 
decomposes according to 


d/iu($) = J n dA * a2 ( a ) / • 

i= 1 


Here, A(A) is the Vandermonde determinant: 


A(A) := det([A^ %) = n( Ai-A ?) 

i>3 


(3.1) 


(3.2) 


and d/iff(Q) is the Haar measure on U(1V). In a partition function, the Vandermonde 
determinant induces a repulsive interaction between eigenvalues: 


Z= j ]JdA; A 2 (A ) J d hh(Q) e -^ s t nAnt ] 

i= 1 

= J fJdAj J d n H (ti) e 


(3.3) 


-/3S[nAnt]+2 E i>7 - log I Ai—Aj I 


In the case of ordinary hermitian matrix models, the action is invariant under T —>• 
nd'Ot. Therefore, the integral over the Haar measure just gives a constant factor. In the 
case of our model (2.10), however, the kinetic term S Un [$] = —^irNtr (d>[a, [o^, <&]]) is not 
invariant under this transformation and hence our first goal is to compute the following 
integral: 

J : = 

As shown in [2], it is possible to rewrite the kinetic term as a sum of traces and multitraces 
of polynomials in under the integral. This multitrace action is then invariant under 
—>• and the integral over H becomes trivial. 

Since our model is invariant under <f> —>• — , our multitrace action will only contain 

terms of even total power of d>. At each order a, there are p(a) terms in Smt of total 


J dp H (ty e /34,rArtr ( nAQt [ a ’[“ t ’ nAnt ]]) (3.4) 
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power a in <5, where p(a) denotes the number of integer partitions of a. We label these 
coefficients by w h ere Tl + ^2 + • • • + is a partition of a: 

Smt[$] =«2 tr (<L 2 ) + opi tr (<f>) 2 + a 4 tr (<f> 4 ) + a 3 ,i tr (<L 3 ) tr (<£)+ 

+ a 2 ,2 tr ($ 2 ) 2 + 02 , 1,1 tr ($ 2 ) tr ($) 2 + 01 , 1 , 1,1 tr ($) 4 + ... . 


Two methods have been developed to compute the coefficients a T1 ,... j7rfc . The first one uses 
group theoretic techniques and was applied in [2] and [3] to compute the partition function 
of scalar field theory on fuzzy <CP N . The second one is a bootstrapping method presented 
in [9], which is more robust and more easily implemented, and we will use this method in 
the following. 

The basic idea behind the bootstrapping method consists in choosing suitable differen¬ 
tial operators D such that 

D e -/3'S'kin[^] — g-^SMT [$] 

D e~P Skin [ <I> ] = e> kin [<£>] e ~ ps ^ [ $ ] and D e -/35MT = O u t [$] e~^ SMT ^ , 


and both O k in and Omt are invariant under <f> —>■ for Q € U(IV). Then the operators 

0 k i n and Omt can be pulled out of the integral and the equation 

D j dfi H (n) = D J e~ fiSMT W (3.7) 


implies O kin = Omt- The left hand side of (3.7) only depends on /3 and N while the 
right hand side depends on /3, N and the coefficients Hence, given a sufficient 

number of differential operators, equations (3.7) will yield enough conditions to fix all the 
coefficients in 5 mt[$], thereby giving the desired rewriting of the action. 

It was found in [9] that the operator dS~ =: dS~ yields conditions that fix more 
than half the unknown coefficients in Smt[^]- More precisely, we have: 

e ~l3S kin m = 16/37riV tr ([&, [a f , $]]) e ~ 0S ^ ] = 0 , (3.8) 


from which it follows that 



e -/3,SMT [$] _ Q 


(3.9) 


This equation holds, in fact, for any scalar field theory on any fuzzy space if the kinetic 
term of the continuum theory vanishes on constant functions and the quantization condition 
<r(l) = 1 is fulfilled. This is the case for our choice of Laplace operator on the fuzzy disc. 

Equation (3.7) for D = yields p(a — 1) conditions on the coefficients of 
which we can use to express all coefficients of the form i in terms of other 

coefficients. In particular, consider the terms in 0 k ; n [<I>] and corresponding to the 

partition tt\ + ... + Ttk-i = a — 1. In terms of coefficients appearing in Smt, Okinfd*] = 
Omt^] amounts to [9] 


a- 


7ri,7T2,...,7r)i,_l,l 


4 ?Z)( (<r(7ri) + 1)£ 


rN 


ri)+l,<T(7r2),...,0-(7T fc _i) + 

+ (<7(^2) + l)®cr( 7 ri),cr( 7 T 2 )-|-l,...,(T( 7 rj,_ 1 ) + • ■ ■ 

T (<7(7Tfc—l) + 1 )®(r(7ri),(T(?r2),...,0'brfc_i)-|-l^ ’ 
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where the sum runs over all permutations of 7 Ti,..., i^k-i and r — 1 is the number of parts 
7 Ti which are 1. Moreover, we define a Wl ,T T2 ,...,n k _ 1 ■= 0 unless tt\ > tt 2 > ... > Ttk—i- I n 
particular, we have 

= (3.11) 

(3.12) 


at second order and 
404 

'Iv 


03,1 = 


- Jj 2 


604 2 a 2,2 


N 


3a 4 02,2 

°id,i,i = -Jp + N2 


at fourth order. Hence, <Smt[$] is determined up to fourth order in 4> by < 22 , 04 and 02 , 2 - 
To fix these, we need to turn to higher order differential operators. 

Unfortunately, none of the higher order differential operators yield functionals dkin)^] 
which are U(lV)-invariant in general. However, simply evaluating the result at 4> = 0 gives 
the desired invariance. To fix 02 , we use D := q§-^q§^- and we readily compute 

d d 


d$ab d'&ba 
d d 

d$ahd<S> 


ba 


which implies 


-/3S kin [«4>] 

— 


$=0 


— 


<4>=0 

a 2 = ~ 

477 N 2 

N + l 


= 8n(3N 2 (N — 1) , 
= -2 (3(N 2 - l)o 2 , 


(3.13) 


(3.14) 


Determining 04 and 02,2 is slightly more involved. We define two differential operators 
D\ and D 2 , 

d d d d , „ d d d d 


D 1 := 


and 


d$ab d^ cd d^ da 
and we solve the pair of simultaneous equations 


D 2 : = 


d^ab < 94 'ba d<5> cd d$dc 


]J. g-^SkinH] _ JJ g—/35 M T [$] 


$=0 


$=0 


* = 1,2 , 


which yields 


04 = 


02,2 = 


8/3it 2 N(12 - 3 N 2 - N 3 ) 

3(1V + l)(lV + 2)(lV + 3) ’ 

8/3tt 2 (36 + 361V - 31V 2 - 111V 3 - 21V 4 ) 


(3.15) 


(3.16) 


(3.17) 


3(!V + l) 2 (!V + 2)(!V + 3) 

The computations for the sixth order coefficients are straightforward but cumbersome 
and we simply quote the result of some computer algebra: 

_ —64/3 2 7r 3 (TV 6 — 15N 5 —5N 4 —123N 3 +388N 2 —30N—120) 

°6 — 3(JV-5)(JV-3)(JV+l)(iV+2)(iV+3)(JV+4)(JV+5) ’ 

_ —64/3 2 7r 3 (2iV 10 +4./V 9 —1037V 8 —811V 7 +1462./V 6 +1610-/V 5 —87837V 4 —58657V 3 +10830iV 2 +2700./V+3600) 
fl4 > 2 — 3(7V—5)(7V-3)(7V-2)7V(7V+1) 2 (7V+2)(7V+3)(7V+4)(7V+5) > 

_ 64/3 2 7r 3 (27V 10 +2TV 9 —987V 8 —1297V 7 + 1634TV 6 + 1384TV 5 —92267V 4 —59057V 3 +139607V 2 +18007V+2400) 
a 3>3 — 3(7V-5)(7V-3)(7V-2)7V(7V+1) 2 (7V+2)(7V+3)(7V+4)(7V+5) ’ 

_ 64/3 2 7t 3 (27V 9 +6TV 8 - 1067V 7 - 1917V 6 +16537V S +30087V 4 -93647V 3 - 157951V 2 +121357V+19020) 
a2 > 2 > 2 ~ 3(7V-5)(7V-3)(7V-2)(7V+l) ii (7V+2)(7V+3)(7V+4)(7V+5) ‘ 

(3.18) 

To keep our computations manageable, however, we limit ourselves to multitrace terms of 
quartic order in 4>. It will turn out that this approximation is sufficient for all our purposes. 
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3.2. Limit of large matrix size 


For computing the partition function of the multitrace action (3.5) plus the potential term, 
we turn to the large N limit in order to apply the saddle point approximation later. Note 
that as usual in quantum field theory, a rescaling of the degrees of freedom needs to be 
accompanied by a rescaling of the fields and the involved coupling constants. This leads 
to a multiscaling limit which we discuss now. 

As N goes to infinity, the discrete set of eigenvalues Ai,..., Ajv of goes over into a 
continuous function A(x) = A(^), 0 < x < 1. Sums over powers of eigenvalues, tr(<J> J ) = 
^N(Aj )- 7 turn into integrals N fg da; A (x) J . Note that each trace yields a factor of N when 
being recast as an integral. 

The maximum total scaling of the terms in the action is ~ N 2 , which is fixed by 
the exponentiated Vandermonde determinant, cf. (3.3). The coefficients of the multitrace 
action scale as follows: 

d 2 ~ N , 04 ^ j3N , (i 2,2 ~ /3 ■ (3.19) 

We denote the scaling of the eigenvalues A, and the couplings /3, r, g by p \, pp, p r and p g , 
respectively: 

P = N p Pfi , X = N px X, r = N pr f , g = N p °g. 

With this notation, we obtain the following inequalities: 

pp + 2 + 2p\ < 2 , 2 pp + 2 + 4 p\ < 2 , 2pp + 2 + 4 p\ < 2 , 

Pp + Pr + 2p\ < 2 , pp + Pg + 4 Px < 2 . 

Each inequality corresponds to a summand in the action. For each summand, we have 
one power of ft outside of the action and further powers from the coefficients, powers of 
N from the various traces and the coefficients as well as scalings of further couplings and 
the eigenvalues. We would like to choose a scaling that saturates as many inequalities 
as possible. The first three inequalities are saturated by pp = —2 p\, and the last two 
equations yield p r = 2 and p g = —2 p\. A convenient choice is therefore 

P\ = Pp = Pg = 0 and p r = 2 . (3.22) 


(3.20) 

(3.21) 


We can now write our model in terms of rescaled quantities. Instead of integrating over 
x , we integrate over A: 

/ dx^ dA p(X) , (3.23) 

Jo Jx 

where p( A) := 4? is the eigenvalue density and X is its support. Introducing the moments 


Ci := dA p(X) X i , 


(3.24) 


we arrive at the action 

/3S = p(rc 2 + gc A - 47t(c 2 - cf) - |/?7r 2 (c 4 - 4c 3 ci + 6 c 2 cf - 3 cf + 2(c 2 - cf) 2 )) + 
- J dXdpp(X) log(|A — p\)p{p) . 


(3.25) 
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4. The phase diagram 

We now come to the computation of the phase diagram of scalar field theory on the fuzzy 
unit disc in the large N limit. We also compare our result to the numerical findings of [11]. 


4-1- Saddle point approximation 
To compute the action (3.25), we rewrite it as 

S\p{X)} = dAp(A)h"(A) - j dAd/ip(A)log(|A - p\)p(p) + £ (^j dXp(X) - 

where the Lagrange multiplier £ was included to fix the normalization of the 
density. The potential reads as 

V (A) = f 3 (aqA ' 4 + Q!3lClA 3 + A~ (c>!2 + Q!21icf + 0^22C2) + A (oqillC 3 + CtllCl)) . ( 4 . 2 ) 


l) , (4-1) 

eigenvalue 


This potential is in fact the general potential for the rewritten action of a quantum scalar 
field theory on an arbitrary fuzzy space 2 , truncated at fourth order in 4>. We recover our 
action (3.25) with the following choice of coefficients: 


an = 47 t , 

16tt 2 /3 
a 22 =- 7 r~ 


aim = 


a 31 = 



«2 = t — 4-7T , 
8tt 2 /3 

«4 = g - 


«211 = 


167r 2 /3 
3 


(4-3) 


The saddle point equation is obtained by varying the action (4.1) with respect to p( A): 


14(A) -2 J dp p(p) log(| A - //|) + £ = 0 , 

where 

V (A) = (5 ^^A 4 + 0:31 (ciA 3 + C 3 A) + CC 2 A'' + o: 2 ii c i A 2 + 

+ 2ci2ii C 2 C 1 A + 2«22C2A^ + 4amic 3 A + 2anCiA 

The key object in finding the solution to (4.4) is the resolvent 


(4.4) 


(4.5) 


H4(A) := J dp ^ , (4.6) 

which is an analytic function on <D\X. A detailed review of the application of the resolvent 
is given e.g. in [13, 14], but for our purposes, the following observations are sufficient. First, 
note that we expect essentially two cases: X can either be given by a single interval or by 
the disjoint union of two intervals. We will refer to these cases as the single cut and the 
double cut solutions. The resolvent’s singular part w(A) contains two roots <5,; in the former 
case and four roots in the latter case. It necessarily satisfies 


cu 2 (A) = M 2 (A) H(X - 5i) = V' 2 (X) - R( A) , (4.7) 

i 


2 assuming that the Laplace operator satisfies Al = 0 
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where M( A) is some polynomial in A and R( A) is a polynomial of one degree less than 
V'(X). The jump over the cut X yields the eigenvalue density according to 

p(A) = “ 2^i ^ (A + k) “ Wl:A “ ’ (4 ' 8) 

which implies that 

p(A) = -^|M(A)| v /((J 2 -A)(A-<y 1 ) , I=[S 1 , 5 2 ] (4.9) 

Ztt 

for the single cut solution or, for the double cut solution, 

p( A) = ^ |M(A)| y/(S A - A)(A - 6 3 )(S 2 - A)(<Ji - A) , X = [ft, d 2 ] U [d 3) <J 4 ] , (4.10) 

Ztt 

where we assumed £4 > £3 > 0 > d 2 > Si. 


4-2. Solutions 


As explained later, we will be interested in two types of solutions. The eigenvalue densities 
of these solutions have support on a single interval and the disjoint union of two intervals, 
respectively. The results of related computations can be found in [15, 16] as well as [2, 3]. 

We will start with the former case and put X = [s — d, s + d\. Equation (4.7) fixes the 
polynomial 

M(A) = mo + mi\ + m 2 A 2 (4.11) 

as well as d. We then use the normalization of p( A), co = 1, and the self-consistency 
conditions for ci, C 2 and C 3 to fix the remaining unknowns. We obtain the eigenvalue 
density 

(4.12) 

(4.13) 


1 


p( A) = —— |mo + miA + m 2 A 2 | y 7 d 2 - (s - A ) 2 

27T 


with 


m0 — /5(2cK2ii<"1 + 3sa3ici + 2a 2 T 2(d 2 + 2s 2 )aq + 4c 2 a 22 ) , 
mi = /3(4 s« 4 + 3ci<a3i) , m 2 = 4/3^4 . 
Additionally, we obtain from (4.7) the equation 

804 s 3 + 6 cia 3 is 2 + (I 2 a 4 (i 2 + 4a 2 + 8 c 2 a 22 + 4c 2 a 2 n) s + 4cian-|- 

+ 3d 2 ci« 3 i + 2 C 3 O 31 + 4cic 2 a 2 n + 8 cfaim = 0 

The normalization condition Co = 1 implies 

\/3d 2 (2a 2 + 2a 211 c 2 + 4a 22 c 2 + 6 a 3 iCis + 3aq (d 2 + 4s 2 )) = 1 , 

and self-consistency conditions for ci, c 2 and C3 read as 

8/3d 2 s (a 2 + 2a 22 c 2 + 3oq (d 2 + 2s 2 )) + 

+ci (/3d 2 ( 8 a 2 ncis + 3a 3 i (d 2 + 8 s 2 )) - 16) = 0 , 


(4.14a) 


(4.14b) 


(4.14c) 
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(4.14d) 


2c 2 (a 22 (3d 2 ( d 2 + 4s 2 ) - 4) + (3d 2 (a 2 ncj ( d 2 + 4s 2 ) + 

+ 6 a 3 icis (d 2 + 2s 2 ) + a 2 ( d 2 + 4s 2 ) + 204 (d 4 + 12d 2 s 2 + 12s 4 ) j = 0 , 

and 

32c3 + (3d 2 ^ — 4 a 2 iic 2 (3d 2 s + 4s 3 ) — 4o 2 (3 d 2 s + 4s 3 ) — 8 a 2 2 C 2 s (3d 2 + 4s 2 ) + 

— 3 a 3 ici (d 4 + 18d 2 s 2 + 16s 4 ) — 12 a 4 S (3d 4 + 14d 2 s 2 + 8 s 4 ) ^ = 0 . 

(4.14e) 

The combined solution of these equations leads to rather involved expressions which do 
not provide any further insight. The special case s = 0, however, does yield manageable 
expressions and it will be of interest to us later on. Here, the eigenvalue distribution is 
symmetric and correspondingly, ci = C 3 = 0. This also implies that a 2 n, « 3 i and aim 
can be put to zero. We obtain the eigenvalue density 


pW 


Vd 2 — A 2 14 — d 2 /3 (d 2 — 4A 2 ) a 4 
2d 2 7T 


where the boundary d is determined by the equation 


(4.15) 


/3 2 a4tt22d 8 + (12/3a:4 + 4/3 o 22) d 4 + 8/3«2d 2 — 16 = 0 . (4-16) 


Note that the above equations can be used to reproduce the results of [3] with the appro¬ 
priate choices of a 2 , a 22 and 04. 

Let us now turn to the symmetric double cut solution with support on the interval 
1 = [— \Js + d, — V s — d] U [Vs — d, Vs + d]. Again, this is a solution with a symmetric 
eigenvalue density, and we thus put ci = C 3 = 0211 = 031 = aim = 0. Following the same 
steps as above, we obtain 


2a 4 /3|A|\/d 2 - (s-A 2 ) 2 

PW = -- • ( 4 - 17 ) 


The normalization condition of the eigenvalue density and the self consistency for c 2 lead 
to 

d 2 /3a.4 = 1 , d 2 /3«4S = c 2 , (4-18) 

which, together with (4.7), yield 


c 2 = s , 


d 


1 

V/3«4 ’ 


a 2 

2 ( 0:4 + 022) 


(4.19) 


4-3. Phase diagram 

Comparing with previous studies of scalar field theories on fuzzy spaces [2, 3] as well as the 
numerical results of [17, 18, 19, 11], we expect three phases: a symmetric single cut solution, 
in which the eigenvalue density has support on a single symmetric interval X = [—d, d], 
a symmetric double cut solution, in which the eigenvalue density has support on two 
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symmetric intervals X = \—\Js + d, — \/s — d] U [y/s — d, Vs + d] and an asymmetric single 
cut solution, in which the eigenvalue density has support on the interval X = [s — d, s + d]. 
These phases are also called disorder phase , non-uniform order phase and uniform order 
phase , respectively. There should be a third order phase transition between the first two 
phases and a second order transition between the last two phases. The former is actually 
the usual phase transition in a hermitian matrix model with quartic even potential, while 
the latter is the analogue of the usual phase transition in two-dimensional scalar quantum 
held theory, cf. [20]. 

We start by considering the existence boundaries of the various phases. We are inter¬ 
ested in the parameter range in which a phase transition can occur, i.e. essentially positive 
g and negative r. Note, however, that for the potential to be bounded from below, we need 
a positive coefficient of A 4 in the potential (4.1). This restricts our parameter space to 


r < 0 


and 



(4.20) 


At the existence boundary of the symmetric single cut solution, the eigenvalue density 
develops a third root at A = 0, signaling the transition to the symmetric double cut regime. 
Putting p(0) = 0, we can solve for d 2 , which, together with the consistency condition Co = 1 
yields an expression for C 2 , which in turn yields the existence boundary 


2 ( a 4 + 0.22) 


(4.21) 


In the case of our model (3.25), the existence boundary is thus at 

2 ( 87 t 2 /3 — g) 

r = — 1 + 47 t . 

(3g - 8 it 2 (3) 

Next, we turn to the existence boundary of the symmetric double cut solution supported 
on the interval X = [—V s + d, —V s — d] U [Vs — d, Vs + d\. Here, one readily finds that 
the existence boundary agrees with that of the symmetric single cut, (4.21). 

Finally, the asymmetric single cut solution only makes sense as long as s — d > 0 and 
p{ A) > 0. This leads to complicated algebraic relations, which cannot be brought into a 
nice analytical expression. Instead, we simply check the validity of our solutions manually, 
whenever required. 

As shown using an approximation in [3], the asymmetric single cut solution as well as 
the symmetric double cut solution exist on overlapping regions of the parameter space. 
To determine the preferred phase, we have to compare the free energy of both solutions 
since a physical system will adopt the phase with the lowest possible free energy. The 
latter is defined as F := — log (Z), where Z is the partition function of our model. In the 
saddle point approximation, we correspondingly have F = /35[p(A)] —/3S'f ree [pfree(A)], where 
5f ree is the free action, truncated at quadratic order in <h and Pf ree (A) is the corresponding 
eigenvalue density. By subtracting the free part, we only let the connected diagrams 
contribute to the free energy. Note that 



PS\p( A)] = J^dX p(X) (V (X) - ±V(A)) - , (4.23) 


11 















as follows from (4.1) and (4.4). The Lagrange multiplier £ is determined by solving (4.4) 
at a suitable point A € X. For example, in the case of the symmetric single cut solution, 
we can choose A = 0 to obtain £ = 2 J x d/i p(/i) log|p| as well as the following expression 
for the free energy: 

F = j|g/3 2 a 2 d 8 + F fi 2 a 2 oi±d & - |/?a 4 d 4 + \f3a 2 d 2 — \ log (-1 d 2 f3a 2 ) — \ (4.24) 


or 


F =lh / 32 ( 3 d ~ 8 tt 2 /3) 2 d 8 + F (r - 4 tt)/3 2 (3 g - 8tt 2 /3) d 6 - ±/3 (3 g - 8 vr 2 /?) d 4 + 


(4.25) 


+ l(r- Ait)/ 3 d 2 - \ log (-irf 2 (f - An) (3) - \ 
in the case of our model (3.25). 

In the cases of the symmetric double cut and asymmetric single cut solutions, we follow 
the choices of [3] and determine £ at A = ±\/i and A = s, respectively. For the symmetric 
double cut, we have 


PS[p(\)] = J d\p(\) 


~ - \ l°g IA — n/s| - 4 log IA + a/sIJ + 

+ \ V$) : 


which evaluates to 


F = — 


Ida 2 


4 ( 0:4 + a 2 2 
and for our model (3.25) reads as 

(3{r — An) 2 


+ 3 log f 


CC 4 


F = — 


4 (ff — 8n 2 (3) ' 4 
For the asymmetric single cut, we have 


+ ilog 


fda? 2 

3 g — 8n 2 (3 
3 (r — 47 t) 2 /3 


PS\p(\)] = j d\p(X) (V(A) - Iy(A) - log |A - s\) + \V{s) , 

which leads to the following lengthy expression: 

Q s 3 «2 2 j6 

F =1§8/3 2 a 2 d 8 + J^s 2 /3 2 ald 6 --f- 3 2 a 2 a 4: d 6 + ^/3 2 c 2 a 4 a 2 nd 6 + 


8 ci 

— i/ 3«4 d 4 + l/3a 2 d 2 — \s 2 f3a\d? + sfdcia^d 2 — 


s^fia^d 2 


+ 


3 ^c 2 a 4 a 2 nd 2 3 s/ 3 cia 4 a 2 nd 2 /3c 4 a 2 n 


4a 22 


4a 22 

„4 / 


+ 


Cl 

s(3c\a 2 


+ i/3c 2 a 2 n(i"+ 


(4.26) 


(4.27) 


(4.28) 


(4.29) 


211 1 


2a 22 2a 22 

3 i 4s 5 /3ct4 


- \ log (- \dj3a 2 ) + 


— ^s z j3a 2 + sf3c\a 2 — 6s 4 /3o4 + 2s 4 /3cia:4 H-41—4 _|_ s/3cian + s/3cfa 2 ii+ 


ci 

i 2 „ 2 /3c?a 2 o:2ii s/ 3 ci 0:20:211 3 s 2 / 3 c?0:40211 3 s 3 /3ci04Q:211 

- ^s z /?c|a2ii H-w-x- i -1-b 


2a 2 2 


2 o 2 2 


a 22 


o 4 0 ,, 3 1 , S I 4sCl c l a 211 , SCI0211 

- pCj^aim + 2spc 1 Q!iiii - 4 + —+ 


«22 

2 s 3 


4 3d 2 3d 2 a 22 d 2 a 22 d 2 3c\d 2 

8 s 4 8 s 3 ci 8 sc 3 a 2 ii 16 s 2 c 2 a 2 ii 8 s 3 ci 02 ii 16s 5 

d 4 3d 4 a 22 d 4 a 22 d 4 « 22 d 4 3cid 4 ’ 


+ 


(4.30) 
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where the variables are subject to the constraints (4.14). 

As consistency checks of our results, we readily verify that the free energies for the 
symmetric single cut and the symmetric double cut solutions agree at the common existence 
boundary, where we expect the usual third order phase transition of hermitian matrix 
models with quartic order potential. Furthermore, the free energy for the asymmetric 
single cut reduces to that of the symmetric single cut if s = 0. Finally, the free energies of 
both the symmetric single cut and the symmetric double cut solutions reduce to those of 
[13, 15, 16] with the correct choice of parameters. 

The phase transition between the asymmetric single cut and the symmetric double cut 
is now found by equating the corresponding free energies. The resulting equations are 
again too involved to be reformulated in analytical form, and we find the corresponding 
solution using simple numerical methods. 

Altogether, we obtain the phase diagram given in figure 1. As discussed above, low 
values of g are forbidden, which is an artifact of our approximation. Also, for very low 
values of both g and |r|, our approximation of the kinetic term becomes unreliable. In the 
remaining parameter space, we have three distinct phases in which the symmetric single 
cut, the symmetric double cut and the asymmetric single cut solutions are appropriate. 
In figure 1, these phases are labeled as I, II and III. The phase transition between I and 
II is the analogue of the usual third order phase transition in hermitian matrix models 
with quartic even potential. The phase transition between II and III is the second order 
phase transition also found in real scalar field theory with quartic even potential on the 
plane. Note that the phase transition from II to III is essentially given by a straight line 
g ~ —4.84r, and computing the phase boundary for high values of — r confirm this feature 
of our plot in figure 1. 



Figure 1: The phase diagram of scalar field theory on the fuzzy disc for f3 = 1. The phases 
I, II and III correspond to the symmetric single cut, symmetric double cut and asymmetric 
single cut solutions. The approximately straight line has asymptotic slope — ^=. 
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Using our formulas, we can now compute the slope of the straight line analytically. Our 
approximation of the kinetic term by multitrace expressions of quartic order in <J> becomes 
better for large values of g and r, and we can restrict ourselves to these. As a second input, 
we find from numerically solving our equations that at the phase boundary where the free 
energies of phase II and III agree, d tends to a fixed value d ~ 4- Both assumptions allow 
us to linearize our equations, which gives us the expressions 

l 

d = _ ~ 0.5250 and g = —^= r ~ —4.84r . (4.31) 

\ 2tt 3V3 

4-4- Comparison to numerical results 

Having computed the phase diagram, let us now compare our results to those of [11]. 

There it was found that neither the coupling g nor the held are renormalized, while r is 

i 

renormalized as r = f N 3. While we also find that only r requires renormalization, our 
factor is different, r = rN 2 . 

The three phases we obtained are also found in the numerical analysis and as stated 
before, the symmetric single cut, the symmetric double cut and the asymmetric single cut 
solutions correspond to the disorder phase, the non-uniform order phase and the uniform 
order phase, respectively. We recover again the fact that the phase transition between the 
non-uniform order and the uniform order phases is given by a straight line. However, the 
best fit of [11] suggested a relation g pc —0.51r, while we found g ~ —4.84r. We are not 
sure about the source of this discrepancy. It might be due to our different renormalization of 
r together with errors introduced by our approximation of the kinetic term by expressions 
of fourth order in or numerical errors for low values of N in [11]. The fact that we 
recover a perfect straight line for larger values of |r| suggests that the approximation of 
the kinetic term is indeed a good one in that parameter range. A last possible cause for 
the quantitative difference of our results to those of [11] might be the different choices of 
Laplace operators. 

Recall that scalar held theory on R 2 with quartic even potential exhibits a phase 
transition with a phase boundary given by a straight line, cf. the numerical results of 
[21]. Scalar held theories on two-dimensional fuzzy spaces necessarily inherit this phase 
transition, and the straight line was found in the numerical studies of [17, 18, 19, 11] on 
the fuzzy sphere and the fuzzy disc. We reproduced this line by analytical methods. 

5. Conclusions 

In this letter, we computed the phase diagram of quantum scalar held theory on the fuzzy 
disc using bootstrapping and matrix model techniques. These methods are analytical but 
employ a perturbative expansion of the kinetic term of the scalar held theory, analogous 
to a high-temperature expansion. The result is a hermitian matrix model with an action 
containing multitrace terms. 
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It had been an initial hope that scalar field theory on the fuzzy disc, which is possibly 
the simplest fuzzy field theory, exhibits some special features allowing a bootstrapping to 
all orders in the field. Unfortunately, this was not the case. 

We proceeded to compute the partition function of the matrix model in the limit of 
large matrix size using a saddle point approximation. We also derived expressions for the 
free energy in various phases. All our results were given in a very general form, applying 
to arbitrary multitrace models with actions of quartic order in the matrix. In particular, 
the free energy for the various phases of quantum scalar field theories on any fuzzy space 
can readily be written down using our equations. 

We computed the shape of the phase diagram for scalar field theory on the fuzzy disc 
with quartic even potential and derived the relevant phase boundaries analytically. We 
find three phases: In the first phase, the disorder phase, the potential has essentially the 
shape of a single well and the expectation value of the magnitude of the field vanishes. 
In the second phase, the non-uniform order phase, the potential is a double well and the 
expectation value of the field distributes symmetrically at the bottom of both wells. In 
the third phase, the uniform order phase, the potential is a deep double well and the 
vacuum expectation value of the field equals one of the minima of the well. Our results 
agree qualitatively with those of [11], but our quantitative results differ. There may be 
various reasons for this discrepancy, amongst which are the different renormalization of 
the couplings as well as errors in the numerical approximations and the different choices of 
Laplace operators. 

The most interesting result of our computations is that we analytically reproduced the 
linear phase boundary between the non-uniform and the uniform order phases. This is a 
feature of both ordinary two-dimensional scalar held theory as well as scalar held theory 
on two-dimensional noncommutative spaces. It was indeed the numerical finding of this 
linear phase boundary for scalar held theory on the fuzzy sphere in [17] which triggered 
the research of [2]. The original aim there was to find the slope of this phase boundary 
analytically. In this letter, we finally achieved this goal and for the case of the fuzzy disk, 
the slope is given in (4.31). 
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